Dynamical Mean-Field Solution for a Model of Metal-Insulator Transitions in 

Moderately Doped Manganites 
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We propose that a specific spatial configuration of lattice sites that energetically favor 3+ or 4 + 
Mn ions in moderately doped manganites constitutes approximately a spatially random two-energy- 
level system. Such an effect results in a mechanism of metal-insulator transitions that appears 
to be different from both the Anderson transition and the Mott-Hubbard transition. Correspond- 
ingly, a disordered Kondo lattice model is put forward, whose dynamical mean-field solution agrees 
reasonably with experiments. 
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The discovery of "colossal" magnetoresistance (CMR) 
has stimulated a renaissance of interest in doped rare- 
earth manganese oxides because of their promising prac- 
tical apphcations and their similarity to the cuprate su- 
perconductor. 1^ Although great efforts have been de- 
voted to this system, there is still no consensus on its 
mechanism so far. In particular, why and how the metal- 
insulator transitions (MIT's) occur in this system for a 
moderate doping range are yet to be clearly understood. 

Various types of phase transitions are displayed in 
the perovskite-type manganites of -Ri-x^xMnOa, with 
R a trivalent rare-earth element and A a divalent alka- 
line earth ion. It is well known that the Mn 3d lev- 
els split in an approximately octahedral ligand into an 
Cg doublet and a t2g triplet; the former is mobile and 
the latter half-filled i2g levels are believed to form local- 
ized spins of S* = 3/2. The parent compound i?Mn03 
is insulating and antifcrromagnetic. Doping of A^^ ions 
creates Mn^^{t2g) holes in a Mn'^+ (tfgeg) matrix; and 
a sufficient number of holes render the mixed-valence 
compound metallic and paramagnetic or ferromagnetic 
depending on the temperature. Besides these doping- 
induced MIT, there is another MIT accompanying the 
ferro- to paramagnetic phase transition as the temper- 
ature is raised in the most interesting doping range of 
0.2 ^ X ^ 0.5. The critical temperature of this transi- 
tion varies with an applied magnetic field, resulting in the 
CMR effect. ||I|,|| Lattice effects [||gl and even structural 
phase transitions p| also interplay, making the enigmatic 
phenomena further intricate. 

The mechanism of the CMR effect is yet controversial. 
It was originally believed to be mediated by the so-called 
"double exchange" since the large ferromagnetic Hund's 
rule coupling tends to align all the d spins as the mo- 
bile eg electrons hop between the Mn ions. |g] Although 
the original model of double exchange in the limits of 
infinite spatial dimensions and classical spins was nicely 
pursued by Furukawa |01 , Millis and coworkers Q argued 
that models involving only double exchange yield results 
such as the magnetic transition temperature Tc and the 
resistivity above Tc that deviate from the experiments 
by orders of magnitude. Accordingly they proposed that 



dynamical Jahn- Teller (JT) effect was crucial. Yet an 
anomalously large electron-phonon coupling of the JT 
origin is used to give a qualitative resemblance to the ex- 
periments PJ. Objections to such vibronic models have 
also been suggested, emphasizing the effects of interfer- 
ence and localization on the electron scattering ||I9 . 

Looking into the electrical, magnetic and structural 
phase transitions, we note the following important facts. 
The doping-induced MIT can appear in both the para- 
and ferromagnetic phases lilt] , and a higher critical dop- 
ing concentration Xc is needed in the former. It is noted 
that a paramagnetic metal has only been observed so far 
in single crystals of Lai_a;Sra;Mn03, whose large tol- 
erance factor (rii;_o/v2?'Mn-o> where r is the averaged 
distance between the two indicated ions.) approaches 1, 
a perfect size match. Special attentions should be paid 
to the sensitivity of the observations to such extrinsic 
factors as preparation conditions [|l2| . 

Accordingly, an essential ingredient is the double ex- 
change leading to the magnetic ordering. This ordering 
in turn facilitates hopping of the eg electrons, resulting 
in a broadening of the electron bandwidth. Another cru- 
cial element is the splitting of the eg bands. In addition 
to the broadening of these bands by the spin ordering, 
which gives rise to a higher Xc for the MIT in the param- 
agnetic phase, they must be broadened by doping as well, 
so that doping can trigger MIT's. The analysis here is 
based on a band-closing MIT, but this is evidenced by 
spectroscopic observations p^ . 

The key point now is the mechanism of the splitting. 
An important fact that is universal in the doped man- 
ganites is the presence of the two valence states, Mn'^+ 
and Mn*"*" , which can even become ordered for appropri- 
ate doping. A salient effect of the doping is to create a 
specific spatial arrangement (random but quenched) of 
v4^+ such that there is a corresponding preferable spa- 
tial distribution of lattice sites that energetically favor 
3+ and 4+Mn ions (we call them 3+ and 4+ sites) @. 
Consequently, that an Cg electron hops from a 3+ sites 
to a, 4+ site to form a Mn^+ will cost an extra energy 
A, which was estimated to fall within the order of mag- 



nitude of the bandwidth even accounting for reasonable 
screening |l5|]. Fluctuations, induced for example by in- 
dividual environment, around the two energies should be 
secondary as the two valence states are prominent. Al- 
though such doping-induced disorder is likely argued to 
trigger a disorder-induced Anderson transition, one usu- 
ally has to invoke the mass-enhancement from the pola- 
ronic effects of size mismatch or of spin disorder in order 
to realize localization nm. We shall show below that 
the main result of the two randomly distributed levels 
is to split an Cg band into two subbands, which can be- 
come overlapped by doping and/or spin ordering. Such 
a physical picture seems to dominate the systems under 
consideration that exhibit conclusively band-closing and 
behave electronically as a linear superposition of its end- 
point compounds |T3]. 

To capture the essential physics, we neglect the fluc- 
tuations around A and include a diagonal disorder with 
only a binary-alloy distribution. As a simplified model, 
we consider only the classical spin limit m. Therefore, 
the Hamiltonian is given by. 
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where mi — {mf, 171^,111^) and |m| = 1, tij 



stands 

for the nearest neighbor hopping integral and J the 
Hund's rule couphng, c\^{cicr) creates (destroys) an elec- 
tron at site i with spin a, and the set of random vari- 
ables Ei assumes an independent identical distribution 
p{e) = x5{e - A/2) + {l- x)S{e + A/2). 

The first two terms of Eq. (P) form simply the Anderson 
model of disorder, while the last two terms represent a 
ferromagnetic Kondo lattice model. Our combination of 
these two models is expected to manifest an interplay 
of disorder, electricity and ferromagnetism. Note that 
a similar but far more complex model that includes in 
addition to the two nonequivalent sites, oxygen orbitals 
as well as strong Coulomb and exchange interactions has 
been used to calculate the band structure pTJ . 

A remarkable reward of the two-level approximation 
is that it renders the model analytically solvable in the 
limit of infinite spatial dimensions d, or, within the dy- 
namical mean- field theory [Q. Although at d ^ 00, this 
method is unable to capture the effects of Anderson lo- 
calization, the results obtained are already nontrivial. In 
this limit, the disordered system reduces to an ensemble 
of self-consistently determined Anderson impurity mod- 
els via proper rescaling of the nearest hopping integral 



tij = tjy^d. The effective bandwidth i, to be set to 1 
below, fixes the energy scale. The on-site Green's func- 
tion G{iujn) is calculated exactly from the effective field 
GoiiuJn) as 



G{iuJn) ^ {{Goiii^n) i-e-Jm-cr) i). 



(2) 



where overbars denote an average over disorder, (• • •)m 
thermal average over m, W and ijj„ = (2n + l)7r//3 = 



{2n + l)TTT are the Matsubara frequencies. On the Bethe 
lattice, the self-consistent condition reads [|l8| 

go{iLLl„)~^ = ILUn + A^ - G{iuJn)/4:, (3) 

ensuring G{iujn) = Y,is.G{'k,iujn) = J D{e)G{e,iuJn)de, 
where D{e) = 1\J\ — e^/vr is the noninteracting density 
of states on the Bethe lattice, and 



G(e, iw„) = T 



1 



iuJn + fJ- — (— T.{iuJn) ' 



(4) 



the exact single- particle Green's function, with the self- 
energy T,{iujn) = GoiioJn)'^ - G{iujn)~^- In Eq.(||), n is 
the chemical potential to be determined by the doping a;. 



l-x = 
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+ 00 



f{uj)A„{uj)duj, 



(5) 



where f{uj) — l/(exp(/3a;) + 1) is the Fermi function, and 
A„{uj) = — ImG'cr(tj+iO+)/7r the single-particle density of 
states (DOS). Moreover, the resistivity is given by p|jlq] 



piT) = polj2jDie)dejA,ie,u;f 



dl{^) 



du) 



duj 



(6) 



where Acr(e, oj) = — ImGo.(e, cu + iO+)/7r and po gives the 
unit of the resistivity. 

For simplicity, only J — > cx) is considered. Defining 
/i = —J + Sfj, and il — uj + 6fi, ||7| and taking the direction 
of the spontaneous magnetization AI to be along the z- 
axis, we finally obtain, (setting cos 6* = ^) 



M= (1-x 



PaMO 



I PAimd^+^l Psmd^, 
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+ 00 



E ln(6on - fci«e ± A) 



(7a) 



, (7b) 



bo = 2n-'-^f^-m^d^-^^ 
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with G„ = AQ - 2bo + 2abi, (ct = ±1), Za^b = 
/expE„ln(&on -fcinC± A)]d^, 6on = bo{iuJn) and 
bin = bi{iuj„). 

For a given splitting A, doping x and temperature T, 
the set of Eqs.(0), can be solved generally by numeri- 
cal methods. Several situations can, however, be fur- 
ther treated analytically. These are T = 0, A = 0, and 



T > Tc, the magnetic transition temperature just below 
which 61 7^ 0. 

For T > Tc, no magnetic ordering appears and so &i — 
0, Pa = Pb and M = 0, and G| = G^ ee g(tj). One 
immediately obtains from Eqs.(0), 



g^ - m-f + (ler^^ 



4A2 + 2)5-8f7 



A{2x- 1)A = 0, 
(8) 



which is identical to Eq. (44) of Ref. 11^ except that 
a different energy scale is chosen. For T — Q, the main 
contribution to Pa and Pb comes from ^ = 1. Conse- 
quently, 61 = 2fi - 60, so G^ = %n- 46o, but Gx = 0, 
i.e., the electrons are completely polarized. One finds 
thus a similar equation to Eq.(@) except with different 
coefficients that lead to a substantial broader bandwidth 
(Fig. [I]). For A = 0, we recover the results of Furukawa 
§. One obtains, for T > T^, g{uo) ^20.- 2i^l/2-D?, 
while at T = 0, G| (w) = 217 - 21^/1-^'^. Thus for mod- 
erate doping, there are only partly filled bands so that 
no insulating behavior exists at all. 

With these known cases as references, it is not difficult, 
albeit tedious, to solve Eqs.(R) self-consistently. Here we 
shall only present the results for the DOS and p{T). We 
choose reasonably the effective bandwidth to be 1.1 2eV 
and po = 10"^ (f7cm)-i. § Thus a splitting A = 0.67 
corresponds to 0.75eV and T = 0.01 to 130.0K. 
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FIG. 1. Density of states for x = 0.2 at different tempera- 
tures. Solid: T = 0.035 (> Tc)\ dashed: T = 0.01 (< Tc); and 
dotted: T = 0. For comparison, the spin down density for 
a; = 0.3 at T = 0.035 (> Tc) (dash-dotted) is given in both 
the main figure and the inset (magnified plot near oj = 0) 



Figure |l| displays the DOS for several values of the 
parameters, manifesting various possible transitions. For 
each orientation of the localized spins, the Cg band splits, 
due to the Hund coupling, into a high-energy antiparallel 
band and a low-energy parallel band; only the latter is 
relevant here, as the separation is of the order of J. Such 
a splitting is similar to the Hubbard coulomb correla- 
tion. The new point here, however, is that for a nonzero 



A, each band splits further into two bands. From the 
spin "down" densities at T = 0.035(> Tc) (identical to 
the "up" ones) of x = 0.2 and x — 0.3 for A = 0.67, it is 
seen that the spectral weight transfers from the lower 
filled electron band to the upper empty hole band as 
X increases. As a result, when the critical doping Xc 
is reached, the two bands overlap, leading to a doping- 
induced MIT. 

At T = 0.01, lower than Tc, the magnetic ordering 
renders the spectral weight for the two spin orientations 
unequivalent (dashed lines). More importantly, it sub- 
stantially broadens the bandwidth, so that the two split- 
ted bands of a; = 0.2 at T > Tc now overlap. Conse- 
quently, an insulator-metal transition occurs accompany- 
ing with the para- to ferromagnetic transition. Moreover, 
as the bands broaden in the ferromagnetic phase, Xc gets 
smaller for the same A. 

On the other hand, as A increases, a larger Xc is re- 
quired. For a sufficiently large A, no doping-induced MIT 
will be present in the paramagnetic phase. This is likely 
the case for systems with a more distorted lattice or a 
smaller tolerance factor. When A becomes so large that 
even for T <C Tc, the two bands are still splitted, a transi- 
tion then occurs from a paramagnetic insulator to a ferro- 
magnetic insulator, or a canted antiferromagnetic insula- 
tor due to the antiferromagnetic superexchange coupling 
between the localized spins. 

The resistivity is shown in Fig. ^ with experimental 
data taken from Ref. pj. The agreements, including the 
positions of the peaks signaling Tc, are remarkable by 
noticing that only a single set of A and po has been 
chosen and several limits have been made in such an ex- 
tremely simplified theory. Note that although for x — 0.3 
and T > Tc the bands have already overlapped, the re- 
sistivity still exhibits somehow an insulating behavior 
{dp/dT < 0). This "pseudogap" behavior has also been 
found in Ref. f^, namely, as soon as a dip, rather than 
a real gap, develops at w = in the spectral function, 
the resistivity starts to rise with decreasing temperature. 
Actually, A and po depend in general on x. So we may 
readily choose other sets of the parameters to produce 
a better agreement with experiments, but this is trivial, 
as Fig. g has already touched the essence. The reason of 
the deviations from experiments is conceivable. As has 
been pointed out above, once the spins are in order and 
hence the bands broaden below Tc, the bonding strength- 
ens, leading to a more compact or less distorted lattice. 
Consequently, the effective bandwidth t increases below 
the transitions, giving rise to the more drastic drops of 
the resistivities observed across Tc. Therefore, a more 
realistic description of the CMR should include the cou- 
pling to the lattice degrees of freedom in order to "dress" 
the "bare" t and also A, so that they can vary with the 
transitions. 
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FIG. 2. Resistivity vs temperatures. Experimental data 
from [2] are shown in the inset. Solid lines: x — 0.2; and 
dashed: x = 0.3. 

In summary, we have put forward a theory that ex- 
plains naturally the doping-induced metal-insulator tran- 
sitions in both the para- and ferromagnetic phases, as 
well as the magnetic phase transitions, with the results 
of doping- and temperature-dependent resistivities being 
in reasonable agreement with experiments. Essential in 
the theory is a random lattice of two-energy levels that 
represent the relative energies of the two species of sites 
that energetically favor 3+ and 4+ Mn ions created by 
doping. So our theory is more appropriate for moder- 
ately doping pfl]. The central physical picture is that 
the Hund's-rule splitted Cg band is further splitted on 
the two levels into a filled electron and an empty hole 
bands, which can become overlapped by doping and/or 
spin ordering. Such a mechanism of electronic MIT ap- 
pears to be different from the two usual classes Q, the 
Mott-Hubbard and the Anderson transitions, since the 
former is due to coulomb correlation, and the latter is 
realized by sweeping the Fermi energy across the mo- 
bility edge. Further, it may be expected to shed light 
to some other multicomponent systems. Although our 
model is extremely simplified, the physical picture for 
the metal-insulator transitions is simple but clear, and 
the theoretical results are in reasonable agreement with 
experiments. The usual Anderson transition should of 
course be investigated. Further extensions of the theory 
are almost transparent. One may study, for example, ef- 
fects of finite J, vibronic coupling, multibands, and the 
Coulomb correlation, to list but a few. Applications to 
other problems are also attracting. 
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